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By combining fermion quantum Monte Carlo (QMC) simulations with diagrammatic theory, we 
have calculated the dielectric screening and the screened electron-electron interaction potential in a 
quasi two-dimensional extended single-band Hubbard model for cuprate superconductors, with and 
without a long-range (1/r) extended Coulomb interaction potential, at and near band-filling i. Our 
results for the insulating dielectric screening in the i-fiUed band case indicate that the Hubbard 
conduction electron system contributes only a minor fraction, Ae = 0.9, of the total observed 
in-plane electronic dielectric constant of the cuprates, e^oW — 4.7. With increasing hole doping 
concentration x, the extended (l/r-) part of the Coulomb interaction potential is rapidly suppressed 
by the metallic screening, due to doping induced charge carriers. Surprisingly, near x ~ 5%, the 
low-frequency component of the screened potential Vs changes sign and becomes attractive at 1st 
neighbor and more extended distances in the CM02-plane. The 1st neighbor screened potential 
reaches maximum attraction strength at dopings a; ~ 13 — 15% and and becomes repulsive again at 
larger dopings, for x ~ 23 — 25%. Similar results are found for the screened potential of the pure 
2D Hubbard model, without extended Coulomb interaction terms, suggesting that the extended 
Coulomb terms do not qualitatively alter the local screening physics of the doped Hubbard electron 
system. When included exactly in the screening calculation, the 1st neighbor extended Coulomb 
term enhances the on-site and 1st neighbor overscreening attraction already present in the pure 
Hubbard model at finite doping. Our results are potentially relevant for the question of the d- 
wave pairing mechanism in the cuprates, since they suggest that the screened extended Coulomb 
interaction potential does not only not suppress a d-wave pairing state but, due to the overscreening, 
could actually increase the d-wave attraction, in the 5 — 25% doping regime. Our results may also 
have implications for the explanation of the isotope effect and its doping dependence in the cuprates. 
At larger dopings, x > 15%, the screened potential becomes attractive even on-site, suggesting 
the possibility that the screened potential could support or enhance s-wave pairing in the high- 
doping regime. We also give a rigorous analytical proof that the screened on-site interaction must 
become attractive near half-filling in the repulsive large-i7-limit of both the pure Hubbard and the 
extended Hubbard single-band model. We present a simple physical interpretation of this seemingly 
paradoxical result in terms of the retardation effects which are inherent in the screening of an 
instantaneous repulsive interaction potential. We also point out that on-site overscreening implies 
singularities in the imaginary-frequency dependence of the irreducible polarization insertion and its 
3-point vertex function. 
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I. INTRODUCTION 

In the cuprate high-Tc superconductors, the strong, on- 
site Hubbard-C/ Coulomb repulsion is believed to prevent 
conventional on-site Cooper pair formation. This long- 
standing theoretical dictum is now supported by exper- 
imental evidence, obtained in several cuprate materiaJi 
for a non-s-wave pairing state, of da;2_y2-symmetry.EJ" 
A d-wave pairing state implies a spatially extended 
Cooper pair wavefunction, involving electrons paired at 
1st neighbor or larger lattice distances. Such a non-s- 
wave state is impervious to the on-site {/-repulsion, but 
it can still be suppressed by the extended part (1st, 2nd, 
... neighbor) of the electron-electron Coulomb repulsion. 
The bare strength of this extended Coulomb potential 



may be quite substantial in the cuprates. For any pro- 
posed microscopic model of the cuprates, it is therefore 
crucial to demonstrate that the model can support a d- 
wave pairing state near i-bandfilling and that such a 
pairing state is robust against extended Coulomb repul- 
sions. 

In the present paper, we combine diagrammatic and 
quantum Monte Carlo (QMC) methods to study both 
the insulating dielectric screening at band-filling ^ and, 
as a function of doping, the screened electron-electron 
interaction potential near the i-filled band limit within 
the framework of a quasi two-dimensional (2D) Hub- 
bard model with three-dimensional (3D) extended l/r 
Coulomb interactions. In the i-filled band case, our re- 
sults suggest that the insulating Hubbard electron system 
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contributes only a minor fraction, Ae = 0.9, to the total 
observed electronic dielectric constant of the cuprates, 
Coo = 4.7, with most of the electronic dielectric screen- 
ing being provided by non-Hubbard "background" elec- 
trons which are not explicitly included in the conven- 
tional single-band Hubbard description of the cuprates. 
In comparing to the observed values for the dielectric con- 
stant, we can thus obtain an estimate for the strength of 
the l/\r\ Coulomb interaction potential of the quasi 2D 
Hubbard electron system. At finite doping, we find that 
screening due to the doping induced "metallic" charge 
fluctuations causes the screened electron-electron inter- 
action potential Vs to become attractive. At low doping 
concentrations x, this "overscreening" effect reverses the 
sign of the spatially extended part of Vs, giving rise to 
a 1st neighbor attraction in the 5 — 25% doping range, 
with a attraction strength maximum near x = 15%. At 
larger x, even the on-site part of Vs turns attractive, for 
X > 15%. We give a rigorous analytical proof for the ex- 
istence of this on-site overscreening effect and show that 
it is intrinsically a doping induced large- 17 effect of Hub- 
bard systems near band- filling i. Possible implications 
of these overscreening effects for superconducting pairing 
instabilities in the 2D Hubbard and quasi-2D extended 
Hubbard model are also discussed. 

The remainder of the paper is organized as follows: In 
Section we describe our model and outline the basic 
approach of combining QMC simulations with diagram- 
matic theory. In Section III, we present our calculation 



of the dielectric constant in the i-filled Hubbard electron 
system and estimate the "bare" strength of the extended 
1/1 r I Coulomb interaction, by comparing to experimen- 
tally observed values in insulating cuprate parent com- 
pounds. In Section IV, we present our results for the on- 
site and 1st neighbor component of the screened Coulomb 
potential in the 2D Hubbard and quasi-2D extended Hub- 
bard model. In Section 0, we present an exact proof for 
the existence of on-site overscreening in the asymptotic 
U oo limit and discuss its implications for the analyti- 
cal structure of the irreducible polarization insertion and 
its underlying vertex function. In Section VI, we spec- 
ulate on possible implications of overscreening for the 
superconducting pairing instabilities of the 2D Hubbard 
and quasi-2D extended Hubbard model. Section VH con- 
tains a brief summary. 



II. MODEL AND COMBINED DIAGRAMMATIC 
QUANTUM MONTE CARLO APPROACH 

Both theoretical (bandstructure) calculations and ex- 
perimental ob- 
servations, specifically the highly anisotropic electronic 
transport properties of the cuprates, suggest that their 
electronic stucture is close to two-dimensiomnal and it 
is a reasonable first approximation to neglect interlayer 
electronic hybridization effects. However, the dielectric 
screening in the undoped insulating parent compounds 



is far less anisotropic, with observed dielectric constants 
of quite similar magnitude parallel and perpendicular to 
the Cu02-layer directions. This suggests that the ex- 
tended 1 /r part of the electronic Coulonb repulsion is by 
no means confined to a single Cu02-layer and that the 
full 3-D interaction potential may have to be taken into 
account. We therefore start from an extended Hubbard 
Hamiltonian of the fornB 
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with ^ creating an electron of spin a ==tj i at Cw-site 
Vj in a three-dimensional (3D) crystal of stacked CuO^ 

layers, rij = X^o- cjo-Cjo- and r^g = Vj —T(. Ht includes 
only an in-plane 1st neighbor hybridization t and the 
chemical potential fj,. Some preliminary results including 
a 2nd neighbor hybridization t' have also been obtained. 
The 3D Coulomb potential 



Vir) = USr,o + 



(2) 



includes the on-site repulsion U and an extended 1/|?'|- 
part, Vc, with a dielectric constant es to account for 
screening by the insulating background not explicitly in- 
cluded in H, that is "non-Hubbard" electrons in lower 
filled bands and, possibly, phonon degrees of freedom. 
On a finite L^^^ x L^^^ x L^^^ lattice with periodic bound- 
ary conditions and primitive unit vectors having 
linear dimension of L^"^ unit cells in the r'^"'-direction 
for a = 1, 2, 3, we define 



mi ,m2 ,m3 



(3) 



a=l 



with the TUa taken over all integers, to ensure proper 
periodicity of V{r). 

The exact screened potential Vs, the exact irreducible 
polarization insertion P and the exact density correlaiion 
function x of the full Hamiltonian are related byO 

Vsiq,iL0) = [1 -Viq)Piq,lu)r'Viq) 

^V{q)-V{q)x{q,^u;)V{q) , (4) 

where V{q) denotes the lattice Fourier sum over V{r) and 
X{q,iu;) = Tdre— '•'■^■MAn,(T)An,(0)) 



-P{q,iu;)[l-V{q)P{q,iu;)]-' 



(5) 



at wavevectors q, Matsubara frequencies iw and tempe 
ature T = l/f3 for lattice size N with Anj = nj — 
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Here, (. . .) and . . . (r) denote, respectively, thermal av- 
eraging and imaginary-time evolution with respect to H . 
Note that Vs and P depend on the choice of representa- 
tion for Hv- Eqs. (Q,^) are based on the diagrammatic 
expansion in the charge representation where the [/-term 
is written as -j , rather than the more familiar spin 

representation U^^nj^jiji. While both are equivalent 
when all diagrams are summed exactly to all orders, the 
former, as we will discuss, may offer some advantages for 
approximate diagram resummations. 

The basic idea of our QMC approach is to calculate 
Vs, via Eq. (^, from the polarization insertion P which, 
in turn, is extracted, via Eq. (|), from QMC results for 
the density correlation function x- I^i order to include 
the full 3D extended Coulomb repulsion into this ap- 
proach, we treat only a certain short-range portion of 
V{r), denoted by Vo{r), exactly, by QMC methods. The 
remaining weaker, but long-range part of V ^ denoted by 
Vf(r) = V{r) — Vo{r), is then handled perturbatively. 

By retaining only in-plane terms in Vo, the QMC sim- 
ulation can be restricted to a single 2D layer. We thus 
calculate, by QMC, the density correlation function 




(Anj(r)„An^(0)o) 



= ~Po{q,luj)[\~Vo{q)Po{qMV (6) 
for the QMC Hamiltonian 

Ho = Hv„ + Ht (7) 

on a single 2D CM02-plane, with the interaction potential 
V in Eq. (|l|) replaced by Vo- Here, (. . .)o and . . . (r)o 
denote, respectively, thermal averaging and imaginary- 
time evolution with respect to Ho and S'o is the exact 
irreducible polarization insertion of H^ which can be 
extracted from the QMC results for Xo via Eq. (||). 

Our essential approximation is to replace the exact P 
of the fuU Hamiltonian H in Eqs. (|,||) by Po, i.e., we set 

P{q, luj) - Po{q, iLu) = [Vo{q) - l/Xo{q, (8) 

in Eqs. and (||). All renormalizations of P due to 
the short-range part of the interaction potential, Vo, 
are thus included exactly, to all orders of Vo- Renor- 
malizations due to the weaker long-range part V{ are 
neglected in P, but approximately included in Vs, via 
Eq. (^, and in X: via Eq. (||), i. e. by setting 
Vs{q,iu;) ^ [1 - V{q)Po{q,i^)]-'V{q) and x(9,H = 
-'Po{q,iuj)[l — V{q)Po{q,iuj)] ^. Note that this approx- 
imation for P becomes exact if we replace H by Ho, as, 
e.g., in our calculations of P and Vs for the pure 2D 
Hubbard model discussed below. 

As a simple cuprate system, we con- 
sider La2-xSrxCu04, with a body-centered tetragonal 
(bet) 3D model crystal structure, in-plane lattice constant 
a = 3.80A, iater-layer spacing d = 6.62AB, t — O.i'beV 
and U = Sot. (Small orthorhombic or larger-unit-cell 



tetragonal distortions from this idealized structure are 
neglected.) Note that the foregoing values for U and t 
are consistent with the values derived by the systematic 
mapping of the low-energy Hilbertspace of the three-band 
onto the single-band Hubbard model.Q The values of U 
and t are also consistent with (and, essentially, deter- 
mined by) the observed values of the antiferromagnetic 
exchange constant J = 0.13eV and of the Mott-Hubbard 
charge excitation gap Amh ^ 1.8 — 2.6eV in the undoped 
La2Cu04 parent compound. □ 

Using standard finite-temperature fermion determi- 
nant QMC methods, Xo is simulated with up to 2 x 10^ 
MC sweeps and typically < 0.5% statistical error on 6 x 6, 
8x8 and 10 x 10 2D square lattices with periodic bound- 
ary conditions at /3 = 1/T up to /3t = 10, with imaginary- 
time step At = (3/ Lr <0.0625t^^ where is the Trot- 
ter number. In most simulations, we use for Ho the pure 
Hubbard model with 

Voir) = USr^o (9) 

In order to explore the effects arising from inclusion of 
extended T^-terms in the polarization insertion Po, we 
have also performed a few simulations for the 1st neigh- 
bor extended Hubbard model, with 

Voir) = U6r,0 + ViY, ^r,s (10) 
seAfi 

where 




is the strength of the bare in-plane 1st neighbor Coulomb 
repulsion and Afi denotes the set of in-plane 1st neighbor 
lattice vectors, s = (±a,0, 0) and (0,±a,0). 

III. DIELECTRIC SCREENING IN THE 
i-FILLED MOTT-HUBBARD INSULATOR 

Our first objective is to calculate the dielectric con- 
stant of the insulator from our model and thereby ob- 
tain an estimate for the background screening es [in 
Eq. (|^)] from the measured long-wavelength external 
field dielectric tensor ^cxU^lq^o in the undoped, insu- 
lating La2CuOi material.D'El Note here that the observed 
eox('^)|g-^o contains contributions from both the conduc- 
tion electron system, explicitly included in our Hubbard 
model description, and from the "background" electron 
and phonon degrees of freedom which are not explicitly 
included in our single-band Hubbard description. Note 
also that in our simplified model, is represented as 
a momentum- and frequency-independent constant, i.e., 
in effect the background is treated as a structureless ho- 
mogeneous dielectric medium, down to atomic length 
scales. We should caution from the outset that such an 
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approach can provide only reasonable order of magnitude 
estimates. 

In a more quantitative modeling approach, based e.g., 
on a shell model representation of the insulating non- 
Hubbard electron and phonon background,u one may in- 
clude more realistically the momentum and frequency 
structure of es. However, as discussed briefly below, the 
"local field" effects in such a model, arising from the mo- 
mentum dependence of es, will not qualitatively alter 
our main conclusions. Also, the frequency dependence of 
the electronic contribution to is likely of lesser impor- 
tance, due to the high charge excitation energy scale of 
the background electrons which is expectecD to be in the 
several eV range, exceeding the Hubbard electron band- 
width 8t. Via the phonon contribution to , on the other 
hand, the low phonon frequency scale fiph ^ O.leV <C 8t, 
enters into the screening of the Coulomb interaction, and, 
along with it, isotopic mass dependence. This may have 
implications for the explanation of the isotope effect ob- 
served in the superconducting properties of the cuprates, 
especially in the underdoped regime, as discussed further 
in Section IV. However, as we will also explain in Sec- 
tion IV , neither the presence of the phonon contribution 



in cb nor its low frequency scale, will substantially al- 
ter the Hubbard electrons' "metallic" screening at finite 
doping which is the central focus of the present paper. 

Within the framework of our simplified background di- 
electric model, we obtain the longitudinal (scalar) dielec- 
tric function e^x fromlij 
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(12) 



at iuj — 0. Here, Vc denotes the 3D unit cell volume and 
as discussed above, x = ^^[1 ^ ^(9)^]^^ from (||) with 
P = Pg. The components of the static dielectric tensor 
fox in the long-wavelength limit can then be extracted 
from 



lim lim 

\q\—>0 L^co 



,{\q\e,iuj = 0) 



(13) 



for arbitrary unit vectors e with L denoting the linear 
lattice size. Note here that it is crucial to take the ther- 
modynamic limit L ^ 00 at finite |g| > before taking 
the long- wavelength limit |g| 0. 

Since the finite-system QMC simulations provide us 
with density correlation data Xo(q, iw) only on a finite, 
discrete g-grid, some special care must be taken to ex- 
tract the foregoing limits from the QMC data. To do so, 
we have developed a simple, physically motivated r-space 
embedding procedure which provides us with a continu- 
ous g-space interpolation and gives a reasonable approx- 
imation of Xo(g,iw) in the thermodynamic limit. Our 
embedding procedure is based on the observation that, in 
the i-filled limit at large U,U ^ 8t, and low temperature 
T, the single-layer density correlation function Xo(r, jw) 
in r-space is actually of quite short range and, already 
on rather small lattices, exhibits very little finite-size de- 
pendence. Physically, this is simply a manifestation of 



the fact that the charge excitation spectrum exhibits a 
large Mott-Hubbard gap, Amh ~ C/ > T. The Mott- 
Hubbard gap not only suppresses all charge correlations; 
it also endows them with a very short charge correlation 
length, ^. In fact, from our QMC data taken on 8 x 8 
lattices at T = O.lt wc estimate ^ to be substantially less 
than the in-plane lattice constant a. 

Based on this observation, we approximate xi^''(r, iw) 
for a large L x L square lattice by the QMC result 

Xo^^'~^\r,iuj), obtained on a smaller lattice of size Lo x 
Lo, for those r- vectors which fall within the (properly 
symmetrized) boundaries of the smaller lattice; for r- 
vectors outside of those small-lattice boundaries, we sim- 
ply set Xo^''(^j*i^) to zero. Specifically, for 2D lattice 
vectors r = (mi a, m2a) on the L x L "host" lattice, with 
integer mi,m2 and |toi|, |m2| < L/2, we set 



ii;(r,L„)x('^^^Hr,zc.) 



(14) 



where, for even 



if \mi\ < Lo/2 and \m2\ < Lo/2 
^ if |mi| = Lo/2 and |m2| < Lo/2 
w{r,Lo)=<^ I if jmij < Lo/2 and |m2| = Lo/2 
if |mi| — Lo/2 and |m2| — Lo/2 
otherwise . 



(15) 



In principle, we could also set the embedding weight fac- 
tors w{r, Lo) to zero for those lattice vectors which fall 
"on the boundary" {i.e. for |mi| = Lo/2 or |to2| = Lo/2), 
without noticeably changing our final results, since even 

for our small QMC lattices, our Xo^^'~^'^ (^j i^) a-re already 
negligibly small on the boundary. With the above choice 

of w, Eq. ([1^), we are ensuring that the large-lattice xi^"* 
provides a "natural" interpolation of the small-lattice 
^(QMC) g-space,in the sense that, after Fourier trans- 



form, identically 



x'o^Hq,^^) 



X 



(QMC) 



{q, iuj) 



(16) 



for those discrete 2D g-points q = {2-k / a){pi/ Lo,P2/ Lo), 
with integer pi,p2, which lie on the discrete g-grid of the 
smaller Lo x Lq lattice. 

An additional complication in extracting eex from the 
QMC data arises from the fact that the QMC simulations 
are performed at finite temperature where the system 
exhibits a small, but non-zero, concentration of electron 
and hole charge carriers due to thermal excitation across 
the Mott-Hubbard gap. In our QMC simulations in the 
grand-canonical ensemble, these thermally excited carri- 
ers manifest themselves by the fact that 



Xo 



(^)((7 = 0,zc^ = 0) >0 



(17) 



i.e. at finite T there is a finite thermal fluctuation in 
the system's total particle number. From QMC simula- 
tions at different temperatures, we have verified that this 
thermal carrier contribution to Xo does indeed exhibit 
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(T) 

the expected activated behavior, i.e., roughly Xo ^ 
exp(— Am/z/T). If included in the calculation of iuj) 
this finite q ~ charge fluctuation would give rise to a 
"quasi metallic" singular contribution to the scalar di- 
electric function which would completely dominate the 
g — !■ limit, with 



f-cx,33 



(24) 



i{q,iuj = 0) 



(T) 



(18) 



Physically this reflects the fact that the thermally ex- 
cited charge carriers are screening out any macroscopic 
Coulomb field with a finite screening length, which is 
simply a manifestiation of the general principle that no 
physical system can be a true insulator at finite temper- 
ature. 

To estimate the (T = 0) insulator contribution to Xo, 
we thus subtract out the thermally activated carrier con- 
tribution and use as our Xo-input into Eq. M) 



(19) 



to calculate Po{q,iuj). With this subtraction, the result- 
ing X from Eq. (^) indeed, by construction, acquires the 
correct "insulating" long-wavelength behavior to give a 
finite Eox for |g| 0. 

Taking the limit L — > oo on xi^iL ^^'^ then expanding 
the resulting xiTns around q = 0, we can write 

X^%{q,^c^ = 0) = e■B■e\q\' + O{\qr) (20) 

for q = \q\e with unit vector e and with the cartesian 
components of the i?-tensor given by 



for a, (3 =1,2,3 
(21) 



„ 1 , ^'x^%icl,^u^ = 0) 
2 igHo oqadqp 



Inserting xiTns('?'*^ ~ 0) ^^ova Eq. ( pO| ) into Eq. (||), 
combining with Eqs. (||,|l^), taking the limit \q\ in 
Eq. ( p^ ) and comparing to Eq. ([l^), we get 



47re^ ~ 

(-bI H B 

V, 



(22) 



With the cartesian coordinate axes chosen along the con- 
ventional tetragonal symmetry directions, only the di- 
agonal components of B and Ccx are non-zero with the 
in-plane component given by 



_ 47re2 

e|| = Ecx,!! — eex,22 = Efi H — — -D| 



(23) 



Since, in our approximation, we do not include any in- 
ter layer charge correlations in Xo^ the out-of-plane ("c- 
axis" ) component of B vanishes and we get 



i.e. the Hubbard electron system does not contribute to 
the c-axis dielectric screening. 

From our QMC results for xi'^^''^*^'' on a io x io = 8 x 8 
lattice, embedded in large LxL lattices with i up to 512, 
we estimate using Eq. (^l|), and from it the Hubbard 
electron system's contribution to the in-plane dielectric 
constant 



II -ej. 



-Bu 



(25) 



Using the pure on-site Hubbard potential, Eq. (||), as 
our interaction potential Vo in the QMC Hamiltonian, 
with the standard Hubbard model parameters stated in 
the previous section, we find Ae = 0.70 from QMC data 
obtained T = 0.333^ and Ae ^ 0.93 from QMC data 
obtained T = Q.lt. Note that this value of Ae is inde- 
pendent of es, since Xo and -B|| depend only on the QMC 
Hamiltonian Ho which is independent of es for the pure 
Hubbard case. 

The slight 25%) T-dependence of the foregoing Ae 
result suggests that the thermal fluctuations, due to ther- 
mally excited carriers, are still noticeably (but not sub- 
stantially) affecting the "insulating" charge correlations 
in the Mott-Hubbard insulator itself, at T = 0.333t. 
This is not entirely surprising since, e.g., the thermal 
carrier subtraction Xo'^'* at that T is about 10 — 20% 
of typical Xo-values at typical non-zero \q\ ^ T^/a, e.g. 

xP I {Xoiq^ioJ = 0))bz = 0.15 where {...)bz denotes the 
Brillouin zone average over q. However, it is reasonable 
to assume that the T = O.U result approximates the 
T = limit to within a percent or better, since the ther- 

(T) 

mal carrier subtraction Xo at T = O.lt is already 5 
orders of magnitude smaller than, say, typical Xo, e-Q-i 
x'-P /{Xo{q,tuj = 0)) BZ = 2 xio-^ 

To explore the effects of including renormalizations due 
to extended Coulomb interactions in the polarization in- 
sertion Po, we have also performed simulations for a ^ - 
filled 1st neighbor extended Hubbard model, Eq. ([iQ), 
with Vi = t = 0.35eV at T = 0.333i. The resulting 
change in Ae, 6Ae = Aely^^t - Ae\vi=o = -0.014 , is 
only about 2% of Ae. Unfortunately, minus problems in 
the QMC simulations prevent us from extending these 
Vi-dependent studies to lower T and/or more realistic 
larger Vi-values, with, say, Vi = 2 — 3 x t. However, 
the foregoing result clearly suggests that the extended 
Coulomb effects on the polarization insertion are entirely 
negligible, at least as far as the long-wavelength limit at 
i -filling is concerned. 

The measured values for the dielectric tensor ecx(t^) of 
undoped La2CuOi in the static limit w — > are eo,|| = 
€\\{ijj = 0) = 30 ± 3 for the in-plane component, and 
eo,± = e_L(a-' = 0) = 25 ± 3 for the c-axis component-Q 
However, these valuesJnclude a large, in fact, dominant 
phonon contribution.Ell 
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The purely electronic contribution to the dielectric 
screening is observed at frequencies uooo ^ 0.5 - leV, 
which are well above the phonon spectrum flp^ < O.leV, 
but at the same time still well below the electronic Mott- 
Hubbard charge gap Amh ^ 1.5 — 2eV. In this fre- 
quency regime, one finds approximately frequency in- 
dependent values of ioo,\\ = ^\\{^oo) — 4.73 for the in- 
plane componenij and eoo,_L = ^^ii-^oc) — 4.56 for the c- 
axis component .a The corresponding anisotropy Aeoc — 
€oo.\\ ~ £oo,± — 0-2 is substantially smaller than our esti- 
mated value Ae = 0.9. However, in our model estimate, 
we are assuming an isotropic background es, whereas, 
in the real material, the background itself could also be 
anisotropic. If we generalize our model to allow for such a 
background anisotropy, we can estimate the background 
in-plane and c-axis components, e^n esx, for our 
model from the experimental uj^-data and our calculated 
Ae. Namely, from the generahzed Eqs. ([23[|2^j25|) , e^n = 
Eooll ~ ^£ — 3.81 and eB± = ej. — 4.57 , corresponding 
to a background anisotropy Ae^ = ~ es± — —0.8 . 

One of the main conclusions from the foregoing anal- 
ysis is that the background degrees of freedom dominate 
the dielectric constant at ^-filling, with the Hubbbard 
"conduction band" electrons contributing only about 
25% to the total electronic dielectric constant, as ob- 
served in the frequency range. The other main con- 
clusion is that both in the background dielectric screen- 
ing and in the full dielectric screening (including Hub- 
bard electrons), the anisotropy is insignificant, compared 
to the actual values of the dielectric tensor components. 
This conclusion holds equally well for the purely elec- 
tronic contribution to the dielectric screening and for the 
phonon contribution, as evidenced by the nearly identi- 
cal values of eo|| ~ 30 and eo_L ~ 25 in the static limit. 
We will therefore in the following continue to work with 
an isotropic background model, with the electronic (uJoo) 
value eb — \/^b\\£b± = 4 and the analogously obtained 
static 0) value es — 27, providing, respectively, 

reasonable lower and upper limits for es- 

From the estimated es=4 (without phonons) or even 
es = 27 (including phonons), one obtains substantial 



1st neighbor repulsion strengths of Vi 



0.95eV = 2.1t in the former and Vi = 0.14eV = 0.4t in 
the latter case. Thus, the extended part of the Coulomb 
potential, Vc{r), appears to be indeed strong enough, 
that it could severely suppress extended (1st neighbor) 
pairing potentials-wliich are commonlyjuioked in both 
phenomenologieaJll3Tl3 and microscopidlao scenarios of 
(i-wave pairing.li^l It is therefore of considerable interest 
to find out how this bare Coulomb interaction potential 
between the Hubbard electrons is modified due to the 
metallic screening generated by the doped Hubbard elec- 
tron system itself. In the next section, we will turn to 
this question. 
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FIG. 1. Screened Coulomb potential Vs(r) = Vs(r, ioj = 0) 
at (a) on site (r = 0) and (b) at in-plane 1st neighbor lat- 
tice vectors r and (c) Eliashberg A-parameters for in-plane 
1st neighbor d^'2_yi and on-site/lst neighbor s-wave pairing, 
all plotted vs. hole doping concentration x = 1 — (jij) at 
fit = t/T = 3.0 and Art = 0.0375 for en = i (full lines) 
and es = 00 (dashed lines). In (a) and (b), results are for 
La2CuOi with bet lattices of sizes 6x6x6, 8x8x8 and 
10 X 10 X 10, with estimated QMC statistical uncertainties in 
Vs of less than O.OleV. In (c), results are based on 8 x 8 x 8 
hct lattice data for Vs, extracted from 8x8 QMC data for Xo- 
Xd has been multiplied x4 for display. 
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IV. QMC RESULTS FOR Vs 

In Fig. |l|(a) and (b), we show results for Vs{r,iuj) at 
iuj — 0, plotted vs. doping x = 1 ~ (rij) for the on-site 
(r = 0) and in-plane 1st neighbor r-vector in the quasi- 
2D extended Hubbard model on the La2Cu04 bet crys- 
tal structure with = 4 and £5=00. Note that the 
latter cb represents just the pure 2D Hubbard model, 
with the extended part of the Coulomb potential, Vc, set 
to zero. The results for Vs were extracted from QMC 
data for the single-plane density correlation function Xo 
of the pure 2D Hubbard model, i.e., with Vo given by 
Eq. (H) , as explained in Section |l[ and with Vs (r, ius) 
obtained from Vs(g,ia;) by Fourier transform over the 
finite-lattice 1st Brillouin zone. Since we are only inter- 
ested in the short-distance behavior of Vs here, we have 
not attempted to embed our 2D Xo-data in larger 2D 
lattices, as was done in the previous Section for the insu- 
lating dielectric function calculations. Rather, all results 
shown are for LqX LqX Lq 3D lattices, using QMC data 
for Xo obtained on corresponding Lq x Lq 2D lattices, 
with Lo = 6, 8, and 10. 

As shown in Fig. |^(b) a small amount of doping sup- 
presses the extended l/|r|-repulsion for r 7^ and causes 
a sign change in the 1st neighbor and (not shown) in the 
2nd and 3rd neighbor screened potential. Thus, Vs(r) at 
short-range in-plane distances r 7^ becomes attractive, 
for X of order 5%, the attraction strength reaches a max- 
imum at a;~ 10 — 14%, and Vs(r) turns repulsive again 
at a; -23 - 28%. 

As shown in Fig. |^(a), the on-site (r = 0) potential, 
while largely unaffected by screening at a; = 0, is also 
rapidly suppressed with increasing x and it also becomes 
attractive at larger doping, near a; = 15%. Over the dop- 
ing range studied, the screened on-site potential varies 
monotonically with x, in contrast to the extended (r 7^ 0) 
part of Vs. 

As a function of es, Vs varies generally monotoni- 
cally, in the directions indicated in Figs. [l](a) and (b). 
If we increase es from 4 to the upper estimated value 
of es = 27, we find that the resulting Vs will be within 
1 — 2% of the eg — CO (pure 2D Hubbard) results shown 
in Figs. 0(a) and (b). Overall, in comparing the quasi- 
2D extended Hubbard model with es = 4 to the pure 2D 
Hubbard model (es = 00), we note that all the forego- 
ing results are qualitatively unaffected by the extended 
Coulomb repulsion. The primary effect of the extended 
Coulomb terms is to move the 1st neighbor Vs slightly 
in the repulsive direction, thereby suppressing the maxi- 
mal 1st neighbor attraction strength in Fig. 0(b) by less 
than 10%. In addition, the extended Coulomb interac- 
tion shifts the "optimal" doping, where the 1st neighbor 
attraction maximum occurs, from x = 10% in the pure 
2D Hubbard model to a; = 13% in the quasi-2D extended 
Hubbard model with = 4. The screened on-site po- 
tential is shifted slightly in the attractive direction, but, 
again, for doping concentrations up to 20%, the effect is 



typically smaller than 10%. 

These results suggest that, despite the substantial 
strength of the bare extended Coulomb potential Vc at 
near-neighbor distances, the local screened potential in 
the Hubbard electron system at finite doping is largely 
unaffected by extended Coulomb interactions. Recall 
from the previous section that the Hubbard electrons 
are largely ineffective in providing dielectric screening in 
the i-filled insulator. By contrast, the results in Fig. |^ 
suggest that the Hubbard electrons' metallic screening is 
very strong, in fact, strong enough to completely over- 
whelm the extended Coulomb repulsion, once a sufficient 
amount of doping and sufficient electronic background 
screening, with x > 5% and es > 4, say, is present. Ad- 
ditional background screening due to phonons will not 
have any significant effect on Vs under these conditions. 

^^From the foregoing discussion, it also becomes clear 
that the presence and frequency dependence of the 
phonon contribution to es will not substantially affect 
the Hubbard electrons' metallic screening at finite dop- 
ing. If one were to include a phononic frequency de- 
pendence in es, with a typical phonon frequency scale 
^^ph ^ 8i, then as a function of iuj, one would find that 
Vs(r, iw) is given roughly by the pure Hubbard result 
(i.e. the dashed lines in Fig. |l|) at frequencies |«a;| <C r^ph 
where the phonons contribute to the background screen- 
ing; and Vs(f, it^) is given roughly by the extended Hub- 
bard result with es =4 {i.e. by the full lines in Fig. |l|) 
at frequencies where the phonons do not con- 

tribute to the background screening. Over most of the 
doping range of interest, e.g., between x ~ 10% and 
X ^ 20%, this additional "phononic" frequency depen- 
dence of Vs, for \iu;\ on the scale of flph, is thus quite 
small, since the e_B = 4.0 extended Hubbard results are 
not very different the es — 00 pure Hubbard results for 
Vs . More importantly, the qualitative features of Vs , such 
as its doping dependence and the overscreening effects at 
finite doping will not be changed by the phononic effects. 
The foregoing argument implicitly assumes that the char- 
acteristic frequency scale of the electronic density corre- 
lations, i.e., of Xo((7, iw), is much higher than flph. From 
earlier QMC studies, it appears that this condition is in- 
deed satisfied, since the spectral weight for density fluc- 
tuations in the Hubbard model extends over a fcequency 
range comparable to the electronic bandwidth.EHI 

By the same arguments, the phononic frequency de- 
pendence of Vs will become increasingly important at 
low doping concentrations where the metallic screening 
of the Hubbard electrons is too weak to suppress the ex- 
tended Coulomb effects. This observation may offer a 
possible explanation for the doping dependence of the 
isotope effect in the superconducting transition tempera- 
ture of the cuprates. The observed isotope exponent^ a 
is generally a decreasing function of doping, in the under- 
doped regime, and becomes negligibly small (compared 
to the classical BCS value of a = i) in the "optimal" 
doping range near x — 15%. This behavior is qualita- 
tively consistent with the notion that, with increasing 
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doping, the phonon effect becomes less and less impor- 
tant in the screened potential. Along with its doping 
dependence, also the overall magnitude of the isotope 
exponent in underdoped cuprates (which can exceed the 
classical BGS valueEJ) has been a long-standing theoret- 
ical puzzle,E3 given the assumption of an electronic pair- 
ing mechanism. It remains to be explored within specific 
pairing models whether the phononic frequency depen- 
dence of the screened Coulomb potential is large enough 
to also explain the observed magnitudes of a in the un- 
derdoped regime. 

Minus sign problems at finite doping unfortunately 
limit our simulations to T > 0.33i. However, at least 
in that temperature regime, we find |Vs| to be increasing 
with decreasing T. This suggests that the overscreening 
effects {i.e. the attraction in Vs < 0) becomes stronger 
than shown in Fig. [| at lower T. 

The presence of both a strong Hubbard-[/ and finite 
doping density x > are crucial for these overscreening 
effects to arise. If we replace the fully Vo-renormalized 
Pq in (Q) by, say, the non- interacting ("RPA") polariza- 
tion bubble Prpa, we also obtain a suppression of Vs{r). 
However, both the on-site (r = 0) and the short-range ex- 
tended part (r 7^ 0) of Vs remain repwlsive in RPA over 
the whole doping range studied here .123 

In the undoped large-?/ system, Vs(r) is reduced rela- 
tive to V{r), by a roughly r-independent factor compara- 
ble to the ratio eex/es, for r 0, i. e. Vs is reduced rel- 
ative to V but retains a repulsive l/|r|-dependence. This 
is expected for the screening of a l/\r\ potential in an 
insulator and confirms the insulating character of the ^- 
filled Hubbard system. The 1 /r long-distance behavior of 
Vs{r) can of course again be traced back to the fact that 
long- wavelength charge fluctuations are suppressed in the 
i-filled insulating state. For x = 0, x(q, iw) ~ |gp 0, 
and hence also P{q, iuj) ^ \q\'^ for |g| ^0. As a con- 
sequence, the screening denominator 1 — V{q)P{q, icu) in 
Eq. remains finite for q ~+ and Vs{q,iLu) inherits 
the singularity of V{q), thus its l/|r|-dependence 

upon Fourier transforming back to r-space. 

By contrast, in the doped system, the screening takes 
on a noticeably metallic character and the r-dependence 
of Vs changes dramatically, even on the finite lattice 
sizes and high temperatures accessible with our QMC ap- 
proach. Already at small doping concentrations, x ^ 5%, 
Vs(r, zw) dies out much faster with |r| than in the un- 
doped case. For example, at "optimal doping" where 
the 1st neighbor Vs is maximally attractive, the in- 
plane 2nd neighbor Vs [at r = (a, a, 0)] is about 4 times 
smaller in magnitude than the in-plane 1st neighbor Vs 
[at r = (a, 0,0)]. 

The metallic screening also strongly suppresses the 
inter-layer Coulomb repulsion, with Vs at the nearest 
inter-layer distance r = (a/2, a/2, c?) being reduced by 
a factor of probably more than 500 relative to the bare 
V at that distance and by a factor of more than 100 rel- 
ative to the in-plane 1st neighbor Vs. Notice here that 



the quasi-2D Hubbard electron system can produce inter- 
layer screening even though our model does not include 
any matrix elements for inter-layer electron transfer. Un- 
like the in-plane components of Vs , its inter-layer compo- 
nents remain repulsive at finite dopings, over the doping 
range explored in Fig. ^ i.e., we do not find any evidence 
for inter-layer overscreening. 

The essential features of Vs are robust against substan- 
tial modifications of the extended Coulomb terms. Note 
here that, in real materials, the simple l/|r|-dependence 
of Vc{rX-.YfSl of course be modified locally by local field 
effects.t2lE3 However, our results do not change by more 
than 20 — 30% if we in- or decrease Vc{r) locally, at 1st, 
2nd, and/or 3rd neighbor distances, by up to 30%, rela- 
tive to Eq. (|^), The latter is a conservative upper limit 
for such local field effects, based cm. the, cuprates' Wan- 
nier orbital and crystal structure. Using different 
layered 3D geometries, such as a Y Ba2Cu^Oj bi-layer 
structure, also does not change the results in Fig. |] by 
more than a few percent. 

We have also carried out QMC simulations with the 
1st neighbor extended Hubbard model, Eq. (p^, at fi- 
nite doping, thereby including extended Coulomb effects 
in our in our polarization insertion P. Surprisingly, we 
find that this actually increases the 1st neighbor attrac- 
tion of Vs and, to a lesser extent, also the on-site attrac- 
tion. For example, if we simulate and analyze according 
to Eqs. (0,|1) an in-plane 1st neighbor Hubbard model, 
with V{r) = Vo{r) given by Eq. ( pi)| ) and a 1st neighbor 
repulsion strength Vi — 0.5i ~ 0.175eV, then the mag- 
nitude of the 1st neighbor attraction in Vs{r, iuj = 0) at 
a near-optimal doping of a; = 14% increases by about 
16% relative to the pure on-site Hubbard model results 
shown in Fig. |l|(c). Notice that this is a comparison of 
two exact results, since we have used the same poten- 
tial V{r) in Eq. (|) as was used as our QMC potential 
Vo(r), both in the pure Hubbard and in the 1st neighbor 
extended Hubbard calulation, i.e. Xo — X ^'^d Pq = P 
exactly. The foregoing result implies that the 1st neigh- 
bor overscreening (in the 5 — 25% doping range) is not 
only robust against extended Coulomb interactions, but, 
in fact, could be enhanced if extended Coulomb interac- 
tions are included exactly in the polarization insertion. 
This result also suggests the interesting possibility that 
extended Coulomb interactions may actually contribute 
constructively to the pairing attraction for extended su- 
perconducting pair wavefunctionsjHJ and_specifically to 



-pairing, as discussed in Section VI 



The effects 

of including extended Coulomb interactions in the polar- 
ization insertion need to be investigated further. Specif- 
ically, the effects of longer-range Coulomb terms (2nd, 
3rd, ...) neighbor need to be explored. 

Only at unphysically large Vi, i.e. unphysically small 
es, does our approach break down, due to charge den- 
sity wave instabilities. These instabilities are signaled by 
l/x('Z, iti^ = 0) — > at some point in g-space. For the pa- 
rameter range explored in the present paper, this occurs 
only for e_B ^ 2.0. 
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V. EXACT PROOF OF ON-SITE 
OVERSCREENING AND ITS IMPLICATIONS 
FOR THE POLARIZATION INSERTION 

In order to see why a sufficiently large U at finite 
doping X = 1 — {jij) must cause an on-site overscreen- 
ing effect, we consider first the pure Hubbard model, 
V{r) = Udr,o, where, from the exact Eq. (^), 

Vs{r,iuj) = U6r,o-U\{r,iLu) , (26) 

with Vs{r,iLL!) and denoting the respective 

Fourier transforms of Vs{q,iuj) and xil^i'^) back into r- 
space. Clearly, the on-site — C/^x-term in Eq. ( p^ ) is at- 
tractive for all iuj, since x(r = 0,z(x'), the autocorrelation 
function of Arij, is always positive. 

At x — 0, i.e. in the i-fiUed insulator, charge fluctu- 
ations are suppressed by the Mott-Hubbard gap. For 
U — > oo, one finds from a large- J7 expansion that 
Xir,iuj = 0)^O{t^/U^), and, in Eq. @, Vsir = OJuj = 
0) = U — 0{t'^/U) > which remains repulsive for large 
U ^ t. The crucial point here is that, at finite |a;| > 0, 
even an infinitely large U can not completely suppress 
the charge fluctuations, since Uj is not conserved at fl- 
nite I a; I, regardless of U. Specifically, the on-site com- 
ponent x(r = 0,zci;) approaches a positive, non-zero limit 
[of 0{x/t), by a simple U = oo scaling argument] for 
U — > oo. Hence, the — [/^X"term will overcome the bare 
[/-term in ( p6| ) and Vs(r = 0, iw) must become attrac- 
tive, i.e. on-site overscreening must occur, for sufficiently 
large U . Also, as a consequence, for sufficiently large U, 
Vs(r = 0, iui = Q) must change sign as a function of x. 

Formally, the suppression of x at ^-filling arises in 
the strong-coupling expansion because the large repul- 
sive U effectively projects out the low-energy Hilberspace 
sector containing only states without doubly occupied 
sites. Upon projection onto the low-energy sector, the on- 
site charge operators Uj become exactly conserved with 
Uj = 1 and Auj = at all sites for U ~* oo, giving x = 
from Eq. (H). At large, but finite U, intersite charge 
transfer processes via the hybridization term Ht, must 
necessarily go through "virtual" intermediate states con- 
taining at least one doubly occupied site. To 2nd order 
in the hybridization Ht , each of the two matrix elements 
exciting from the low-energy sector into this high-energy 
sector and back is of order t/U while the inverse energy 
denominator associated with the virtual excitation is of 
order 1/U, resulting in x ~ t^//7'^. At finite doping, 
on the other hand, the projected nj are not conserved. 
Even for U — oo, the doping induced holes can still move 
through the lattice, via the Ht-term, since the intersite 
charge transfer can proceed without having to go through 
high-energy intermediate states involving double occu- 
pancy. Hence, x approaches a non-zero limit, Xoo 7^ 0, 
at finite doping concentration for U 00. 

Note that the foregoing large-C/ argument constitutes 
an exact analytical proof of on-site overscreening in the 



asymptotic limit C/ — > cx). The proof holds both on fi- 
nite and infinite lattices in any spatial dimension and 
it immediately generalizes to the full extended Hubbard 
model, Eqs. (|l|,y) since the — C/^x-term remains the dom- 
inant screening contribution for U/t — *■ 00 at finite a; > 0, 
even in the presence of a finite extended interaction term 
Vcj^O. Our QMC results in Fig. |l|a not only confirm these 
exact large-C/ results; they also show that the large-C/ 
scenario is realized in the physically relevant parameter 
regime^ C/~8- 12t. 

The occurence of on-site overscreening, or close prox- 
imity to it, is accompanied by profound effects in the 
polarization insertion P{q,iuj), namely, by g-dependent 
singularities of P on the (analytically continued) iw-axis. 
To see this, recall that Vs(r = Q,iuj) is just the Bril- 
louin zone average of Vs(g, iw). Hence, in order to get 
Vs(r = 0,iijj) < 0, there must be some region in q- 
space, for which Vs{q,io-!) < 0. On the other hand, since 
x{q, ii^) ~ Oil/uj"^) for \iuj\ 00, it follows from Eq. (^) 
that Vs[q,iijj) — + V{q) > for \iuj\ — > 00, i.e., in the 
high frequency limit the screened potential approaches 
the repulsive bare potential. Hence, for those q for which 
Vs{q,iuj) < at some (low) frequency icu, there must 
exist at least one iuj^^\q), on the imaginary frequency 
axis, such that Vs{q, iuj), analytically continued onto the 
continuous icu-axis, goes through zero, 

lim Vsiq,iuj) = 0. (27) 

iuj — ncj^^^ (q) 

By Eq. (||), this implies 

lim l/xiq,tuj)^V{q) (28) 

iuj — >iLJ^^' (q) 

and, by Eq. (^ which is equivalent to \lP{q,iuj) = 
V{q) ~ l/x{q,i^), this implies 

lim 1/P{q,iuj) = . (29) 

In other words, \/ P{q,iuj) changes sign and P{q,iuj) 
must be singular at iuj^'^^q). 

Note that this is clearly a strong correlation effect. 
Weak couphng approximations to P{q, iuo), notably RPA, 
do not give such a singularity in P{q,iu!). In RPA, 
1/P{q, iuj) < at all q and iuu. It is therefore not surpris- 
ing that RPA cannot reproduce the on-site overscreening 
effect. Note also that this singularity in P does not im- 
ply any singularities (i.e. instabilities) in physically ob- 
servable quantities, such as x(9i*^)- In fact, xl?, *'j-') is 
perfectly regular at iu)^^''{q), as implied by Eq. ( p8| ) and 
V{q)>0. 

However, the singularity of P{q, iuj) does imply sin- 
gularities of certain vertex functions to which both 
P and the single-particle self-energy E are diagram- 
matically related. For example, both P and E can 
be expressed in terms of an appropriately defined 3- 
point vertex function A(fc, iv; q, ico), with entering boson 
momentum- energy (q,iuj), entering fermion momentum- 
energy (k,iv) (where ii/ is an odd Matsubara frequency) 



9 



and exiting fermion momentum-energy (k + q^iiy + iu), 
such thatcl 



P[q,iu) 
T 
N 



G{k, w)G{{k + q,iv + iuj)A{k, iv; q, ioj) (30) 



k.iu 



and 



E(fc, iv) 

= SH(fc) 

T 



(]v) X!^s(g,«^^)G(fc-l-g,w + it^)A(fc,w;g,it^) (31) 



where G{k, iv) is the single-particle Green's function and 
Eh denotes the self-energy contribution from the Har- 
tee diagram. Note that the singularity of P implies 
an analogous singularity in the iw-dependence of A for 
ibj iuj'^'^^q). However, in the foregoing expression for 
S, which is, in principle, a physically observable quan- 
tity, the singularity of A is cancelled by the vanishing 
of the other factor in the summand, Vs{q,ii-u) for 
iuj — > itt'*-*-'((7). The formal implications of this singular- 
ity in P and A need to be further investigated. 

As discussed in Section [V, RPA results show that, for 
near (1st, 2nd, ...) neighbor r's, x(r, iw = 0) is negative 
at small U, whereas QMC results suggest that the near- 
neighbor x(r, = 0) becomes positive at finite doping 
when U exceeds a doping dependent threshold of order 
several t. ^^From Eq. ( p6[) one sees that this positive 
near neighbor x(r, iw = 0) gives rise to the near neighbor 
attraction in Vs at large U and finite x, as displayed in 
Fig. I^b. Hence, the QMC and RPA results taken together 
suggest that near-neighbor overscreening is fundamen- 
tally also a large-?/ effect, just like on-site overscreening. 
However, note here that x(r, itu) for r ^ is not an au- 
tocorrelation function and its sign is allowed to be either 
positive or negative, depending on the model parame- 
ters. Because of this, it does not seem to be possible to 
generalize the above large-C/ overscreening proof to show 
the existence of near- (1st 2nd, ...) neighbor overscreen- 
ing in Vs{r,iuj) by analytical means. Also, unlike the 
on-site overscreening case, the existence of near-neighbor 
overscreening does not necessarily imply the existence of 
singularities in P{q, ilo). 



VI. D- AND S-WAVE PAIRING STRENGTHS 

Given the attractive nature of the screened potential at 
finite doping, it is tempting to ask whether this attraction 
could give rise to superconductivity and, if so, of what 
pairing symmetry. A potential advantage of our diagram- 
matic expansion in the charge representation is the large 
reduction of the overall strength of Vs in the 10 — 20% 
doping range, compared to the bare Hubbard-[/. This 



suggests the possibility of carrying out controlled, self- 
consistent weak-coupling expansions in which the fully 
screened Vs, rather than the bare V or U , serves as the 
small parameter. j-Such an expansion can be formulated 
diagrammaticallyD by retaining only skeleton diagrams 
in which none of the interaction lines contain any po- 
larization insertions and each interaction line represents 
a Vs((Z, iw). Superconducting instabilities can then be 
studied in terms of such a perturbative approximation to 
the irreducible particle-particle vertex, expanded to 1st 
order in Vs. 

As a first step in that direction, we have explored 
possible Vs-induced or -enhanced superconducting pair- 
ing instabilities, pusing the standard Eliashberg-McMillan 
(EM) approachnJ A convenient measure of the pairing 
strength of Vs are the dimensionless EM A-parameters, 
defined in terms of the Fermi surface "expectation val- 
ues" of Vs(fc — fc'jiw = 0) for relevant Cooper pair trial 
wavefunctions r){k) in electron momentum (fc-) space, as 
described, e.g., in Refs. |l3| and 

In Fig. 1(c), we show the EM parameters A^, for on- 
site s-wawe (and, identically, for in-plane 1st neighbor 
s-wave),c3 and A^, for in-plane 1st neighbor d^2_y2 pair- 
ing, with respective pair wavefunctions ris{k) = 1 and 
Vdik) — cos(afca;) — cos{aky). To carry out the required 
Fermi surface integrals, our 3D Vs{q, iw) was interpolated 
from the finite 8x8x8 lattice g-grid onto a 200x200x200 
q-grid, using the 3D version of the g-interpolation scheme 
described in Section HI for the 2D g-interpolation of Xo- 
Applying this interpolation scheme to Vs((j, icj) is jus- 
tified here, analogous to the above Xo-interpolation, by 
the fact that, at finite doping, Vs{r,iuj) is of very short- 
range in r-space, due to the "metallic" character of the 
screening. 

At low doping, the dominant attractive (A > 0) channel 
is d^2_y2 with Ad reaching a maximum of ~0.15 — 0.17 
near x ~ 10 — 14%. As is repulsive at low doping, but 
becomes strongly attractive at larger doping x 15%. 
Thus, as expected on symmetry grounds, A^ and A^ re- 
fiect the doping dependence of the on-site and 1st neigh- 
bor attraction Vs shown in Fig. |l|(a) and (b), respectively. 
We note in passing that the doping dependence of A^ is 
reminiscent of the observed dopiru| dependence of the su- 
perconducting Tc in the cuprates.cl The A- values for near- 
neighbor pair wavefunctions of other symmetries (p, d^y, 
g) are small compared to A^ and A^ and for that reason 
not further discussed here. 

The spectraL weight of xi^i^^) extends up to values 



f^x 8 - 10t.E°l In the EM analysisP this "boson" en- 
ergy scale, together with A, determines the superconduct- 
ing Tc, roughly as ~ fi^ exp(— 1/A). Because of the 
large fi^-scale, it may be possible, at least within the 
EM approximation, to achieve high T^'s even for moder- 
ate coupling values A < 1 . 

We should caution here that the foregoing QMC results 
can of course only suggest the general, qualitative trends 
in Vs and A parameters, due to the high temperatures 
and small lattice size limitations inherent in the QMC 
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approach. Also, our EM approach is based formally on a 
perturbative expansion of the exact irreducible particle- 
particle vertex to 1st order in Vs- This approximation 
should be relied upon, if at all, only close to the "cross- 
over" xc — 15% where Vs(^ — 0,iuj = 0) = 0, so that 
|Vs(9:*'^)| *C 8t, at least for low frequencies \ioj\ <C 8t. 
i.e. roughly in the 10 — 20% doping range. At large 
over- or underdoping (x ^ 20% or a; ^ 10%), where 
I As I ^ 0{1), corrections of higher order in Vs are likely 
to contribute strongly to the self-energy and to the irre- 
ducible particle-particle vertex, thereby causing the EM 
approximation, with Vs as the effective pairing potential, 
to break down. The increasingly attractive As-values at 
X ^ 25% doping, in Fig. Ic, do therefore not necessarily 
imply increasingly strong s-wave pairing tendencies. The 
effects of vertex corrections, of higher order in Vs, need 
to be further explored, both in the over- and underdoped 
regimes. 

It is also important to realize that the strong reduc- 
tion in overall screened potential strength |Vs| is only 
a necessary, not a sufficicient condition for the appli- 
cability of the EM approach. Note in particular, that 
Vs is "weak" (and, where applicable, attractive) only at 
low frequencies. At large frequencies, with \iuj\ of the 
order of the bandwidth 8t, the charge correlation func- 
tion begins to die out, with xilji^) and, by 
Eq. (^), Vs('?, *'^) V{q) for \iuj\ oo. In other words, 
at high frequencies, the screened potential recovers the 
full strength of the bare potential. Physically, this sim- 
ply reflects the fact that at high frequencies, ^ 8t, 
the conduction electron system is "too slow" to provide 
a screening response to very rapidly varying fields. Since 
the characteristic frequency scale of these charge fluc=. 
tuations is quite high, of the order of the bandwidth,L2l 
the conventional underpinnings of the EM approxima- 
tion, such as the Migdal theorem, are not necessarily sat- 
isfied. Thus, even in the "near-optimal" doping regime 
where the low-frequency Vs is weak, it needs to be re- 
examined whether corrections to the irreducible particle- 
particle vertex of higher order in Vs are indeed small. 

With these caveats in mind, we should compare 
our results for the pairing strengths in the Hub- 
bard model to earlierj-studies of the superconduct- 
ing pairing cprmlatiaiiscZrE3 and of the effective pair- 
ing potentiaMllZlEl'y in the Hubbard-Jiiodel. Early 
studies of pair correlation function£l£p-| and pair 
bound state symmetries on small lattice^^-^ using for- 
mally exact QMG simulation methodscZrEj^LJ and exact 
diagonalization,E3 indicated a tenjdeacy towards dx'^-y'i- 
pairing. In finite-T simulationsj£3~E3 this was suggested 
by an increase in d^2_y2 pairing correlations with de- 
creasing temperature. However, due to the QMC minus- 
sign problem, these exact QMC studies were limited to 
small lattice sizes and rather high temperatures (in phys- 
ical units about an order of magnitude higher than the 
observed scale in the cuprates). It may therefore be 
difficult to draw meaningful conclusions about the low- 
temperature long-range pairing correlations of the model 



from the finite-system and/or finite-T exact QMC data 
for pair correlation functions. 

Also, by contrast, more recent studies of pairing corre- 
lations in the groundstate of the Hubbard model, based 
on the "constrained path Monte Carlo" (CPMC) ap- 
proach, did not provide any evidence for long-range pait^ 
ing correlations of either dx2-y2- or s-wave symmetry.E2l 
However, the CPMC approach is based on a (in prin- 
ciple uncontrolled) variational approximation which, al- 
beit remarkably successful in reproducing exactly known 
groundstate energies, may not necessarily reproduce the 
exact pairing correlations. Thus, existing QMC studies of 
pair correlation functions in the Hubbard model are, at 
best, inconclusive as far as (i2;2_j^2 -pairing is concerned 
and, with all the above-cited limitations, they do not 
provide evidence for s-wave pairing in the near-^-filled 
Hubbard model. 

In addition to the above-cited "technical" limitations, 
there is a further, physical reason why especially long- 
range s-wave pairing correlations, if existent in the Hub- 
bard model, may escape detection in QMC simulations of 
pair correlation functions. To understand this, note that 
retardation plays a central role in the physical origin of 
the overscreening attraction in our Vs. For example, for 
the pure Hubbard model, Eq. (^6|), Fourier transformed 
to the (imaginary) time (r) domain, becomes 

Vs{r = 0,T) = U5(t) - U\ir^O,T) . (32) 

A pair of electrons, e and e, say, is subject to the instanta- 
neous, bare repulsive U interaction term only if they both 
occupy the same site j at the same time. By contrast, the 
attractive screening term, — C/^Xj is retarded, i.e., phys- 
ically speaking, the charge polarization caused by e at 
site j can still be felt by e after a time lag 9 > 0, when e 
has already left j. Thus, the on-site overscreening arises 
from processes wherein the two electrons interact on the 
same site, but at different times, via their local charge 
polarizations, thereby evading their bare, instantaneous 
repulsion. The time scale for this retarded interaction 
is governed by the frequency spectrum of the dynamical 
density correlation x^ 6-: for the relevant parametjeu 
range it is of order of the inverse bandwidth {8t)^^EB 
Pair correlations induced by such a potential, with re- 
tarded attractive and instantaneous repulsive contribu- 
tions, should be most easily detected via time-delayed 
pair creation order parameters of the general form 

A, = c]je)cl^ (33) 

where c] ^{6) = e^^cj „e^^^ . In other words, Ag creates 
the second electron e with time lag 9 relative to the first 
electron e. For the on-site s-wave (j — j) case, one should 
thus consider time-delayed pair correlation functions of 
the general form 

Cir,-r,,,T,e,0') = {T[c,^ir + e)c,^iT) c],^(0')4^(O) ]) 

(34) 
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where T[...] impHes fermion time ordering. The strongest 
signal for pairing correlations should then be detected for 
typical values of 9 and 9' which are comparable to the re- 
tardation time scale of the attractive part of the poten- 
tial, i.e., in the case of our Vs(r, r), for 9, 9' ^ {8t)^^. By 
contrast, the existing QMC studies of pairing correlation 
functions in the Hubbard model have so far considered 
only simultaneous pair correlations C{rj — rji ,t,9,9') 
with 9 = 9' — 0. In such a simutaneous pair correla- 
tion function, the detectable signal for s-wave pairing is 
expected to be very weak (and, in the case of pairing 
with odd frequency parity, identically zero) if the under- 
lying pairing potential happens to have a time-dependent 
structure of the sort described above for Vs(r, r), i.e. 
a strong instantaneous repulsion superimposed on a re- 
tarded attraction. We note that for such a temporal 
structure of the pairing potential the possibility ofpa 
triplet s-wave pairing state with odd frequency parityE3 
can not be ruled out. 

In using QMC simulations to find superconducting in- 
stabilities via the pair correlation function approach one 
thus faces at least three major difficulties: (i) long-range 
pair correlations develop only at low temperatures, of or- 
der Tc, (ii) it requires large lattice sizes to detect such 
long-range correlations in finite systems (a problem com- 
pounded at finite T in two dimensions by the expected 
algebraic decay of the correlations); and, (iii) depending 
on the structure of the superconducting order parameter, 
the simulated pair correlation function may provide only 
an undetectably weak signal if one chooses the "wrong" 
pair creation operators (and one doesn't know a priori 
which operators are the "right" ones). A potentially 
promising route to circumvent some of these difficulties is 
to focus instead on the underlying exact effective pairing 
potential, i.e. in precise diagrammatic ter pi s ,jpii jthe ir- 
reducible particle-particle vertex function.taOO'ES The 
basic idea here is that, in contrast to long-range pair 
correlations, the effective pairing potential may be (i) of 
short range and may (ii) develop strong signatures of the 
pairing attraction already at high temperatures. Also, 
(iii) by solving for the dominant pairing eigenvalues and 
eigenf unctions of the the corresponding particle-particle 
ladder equation, the particle-particle vertex approach au- 
tomatically reveals the symmetry and space-time struc- 
ture of the dominant superconducting order parameter. 

The existing QMC results for the exact irreducible 
particle-particle vertex of the 2D Hubbard model near i- 
fiUing suggest that there is indeed a noticeable 1st neigh- 
bor pairing attraction and that a pairing eigenfunction 
of spin-singlet d,j.2-y2 symmetry becomes dominant {i.e., 
develops the largest .pairing eigenvalue) as the temper- 
ature is lowered.EStJ In addition, sub-dominant pairing 
eigenvalues of comparable magnitude with eigenfunctions 
of odd-frequency triplet s-|-apd odd-frequency singlet p- 
wave symmetry are foundJlj The effective "on-site" in- 
teraction extracted from the irceducible particle-particle 
vertex is found to be repulsive .liZI 

We should emphasize here that, aside from some pre- 



liminary studies at ^-fiUingJSS all of the foregoing exact 
irreducible particle-particle vertex results were obtained 
at a doping concentration of x = 13% and, at that con- 
centration, they are not inconsistent with our results for 
the pairing A-parameters shown in Fig. ^c). At a; = 13%, 
our dj.2_y2-A is indeed the dominantly attractive one; the 
s-wave-A is repulsive at x — 13% and becomes attractive 
only for x > 15%. 

Also, unfortunately, the reported results for the real- 
space interaction strengths extracted from the irreducible 
particle-particle vertexliZI were reported for smaller Hub- 
bard U,U — At, whereas a physically more realistic larger 
value ofU — 8t was used in our calculations of Vs(r, iui). 
Recall here that, as discussed above, the on-site over- 
screening is a large-f7 effect and U — At may simply be 
too small to produce on-site overscreening. Interestingly, 
from a comparison of the pairing eigenvalues extracted 
from the irreducible particle-particle vertex at C/ = 4i 
and at U = 8t,E3 it appears that the (odd frequency 
tripet) s-wave eigenvalues are increasing with U. This 
observation is not inconsistent with the notion that one 
moves towards on-site overscreening with increasing U. 

It would therefore be of considerable interest to check 
whether the irreducible particle-particle vertex develops 
an on-site attraction and stronger s-wave pairing ten- 
dencies at larger doping, x > 15%, and larger Hubbard 
repulsion U, e.g., for U — 8t. This would provide a very 
stringent test as to whether our screened potential Vs in 
the charge representation indeed provides a reasonable 
approximation to the exact irreducible particle-particle 
vertex or whether particle-particle vertex corrections of 
higher order in Vs are important. 

In analyzing the irreducible particle-particle vertex, 
it was also found that its momentum and real-space 
structure, and the T-dependence of that momentum 
and real-space structure, closely resembles that of an 
appropriately defined spin fiuctuation mediated pairing 
potential of the Hubbard model, suggesting that spin 
fiuctuatons ar e i n deed responsible for the d-wave pair- 
ing attractionOiy It is possible that the positive near- 
neighbor charge correlations [i.e. x(r, = 0) > at 
near-neighbor r- vectors] , responsible for the d-wave pair- 
ing attraction in our approach, are closely (and per- 
haps causally) related to short-range antiferromagnetic 
spin correlations in Hubbard systems near ^-filling. The 
charge representation approach developed here may thus 
provide a description of the physics in near-i-filled Hub- 
bard systems which is comjplpmetitary to that of a spin 
fiuctuation-based approach J3l2a^EZl The overscreening of 
the on-site potential, and hence the possibility of s-wave 
pairing in the Hubbard model, is one aspect of this prob- 
lem which may be (with all the above-stated caveats !) 
"obvious" in the former, but difficult to reproduce in the 
latter approach. The relationship between our charge 
representation formulation -and the spin fluctuation pic- 
ture of the Hubbard modcHO'EZl needs to be further in- 
vestigated. 
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VII. SUMMARY 

In summary, by a combination of diagrammatic and 

quantum Monte Carlo techniques and by exact analyt- 
ical approaches, we have studied the dielectric screen- 
ing and the screening of electron-electron Coulomb re- 
pulsions in the charge representation of the 2D Hubbard 
and of the quasi-2D extended Hubbard models relevant 
to the cuprate superconductors. 

Applying a combination of QMC simulations and di- 
agrammatic techniques, we have calculated the contri- 
bution of the quasi-2D extended Hubbard electron sys- 
tem at ^-filling to the long-wavelength in-plane dielec- 
tric screening in the undoped insulating cuprate parent 
compoimds. Combining this result with experimental 
data for the observed long-wavelength dielectric tensor 
in the La2CuOi, we have then obtained an estimate for 
the strength of the dielectric screening due the insulat- 
ing "background" of non-Hubbard electrons and phonons 
and, thereby, for the effective strength of the extended 
l/|r| Coulomb interaction potential in the Hubbard con- 
duction electron system. Our results imply that no more 
than 25% of the observed electronic in-plane dieclectric 
constant (cqo) arises from the ^-filled Hubbard conduc- 
tion electron system. 

Using again a combination of diagrammatic and QMC 
techniques, we have then studied the effective, screened 
electron-electron interaction potential Vs, as a function of 
doping X, in the both the 2D Hubbard and quasi 2D ex- 
tended Hubbard model. We find that finite doping gives 
rise to overscreening effects which cause Vs to become at- 
tractive at finite doping concentrations. Specifically, we 
find that the low-frequency extended (1st, 2nd... neigh- 
bor) part of Vs, while repulsive and rapidly suppressed 
with doping at small x, becomes attractive for x > 5%. 
Near a; ~ 10 — 15%, the extended part of Vs exhibits 
maximum attraction strength and, for .t 20 — 25%, it 
becomes repulsive again and increases with x. At larger 
doping, x> 15%, even the on-site part of Vs changes sign 
and becomes attractive, with a monotonic ^-dependence 
between a; = and x = 30%. 

Both spatially extended (1st, 2nd, ... neighbor) and 
on-site overscreening effects are robust against the pres- 
ence of and moderate changes in the 3D extended l/|r| 
Coulomb repulsions, and largely independent of the ar- 
rangement of Cu02-layers in the 3D crystal structure. 
Generally, at doping concentrations larger than 5 — 10%, 
our results for the screened interaction potential of the 
extended Hubbard model, with extended l/|r| Coulomb 
interactions of realistic strength, differ from those of the 
pure 2D Hubbard model (with only on-site U repulsion) 
by no more than 10 — 20%. This suggests that the ex- 
tended Coulomb effects do not qualitatively alter the the 
short-range charge correlations and screening in the Hub- 
bard electron system at finite doping near ^-filling. In 
particular the overscreening effects which we find at dop- 
ing concentrations x > 5% are intrisic to the pure 2D 



Hubbard model. These overscreening effects arc modi- 
fied, but neither caused nor destroyed by the extended 
Coulomb terms. When included exactly in the screening 
calculation, a repulsive 1st neighbor extended Coulomb 
term enhances the 1st neighbor overscreening in Vs at fi- 
nite doping, suggesting the possibility that repulsive ex- 
tended Coulomb terms could actually contribute to an 
extended pairing attraction. 

Independent of the QMC results, we have presented an 
exact analytical proof that, in the large- [/ limit, doping 
induced charge fluctuations must give rise to an "over- 
screening" effect which causes the exact screened on-site 
Coulomb potential Vs in the charge representation to be- 
come attractive. We have shown that this is indeed an 
intrinsic local property of both large-?/ Hubbard and ex- 
tended Hubbard models at finite doping which exists in- 
dependent of dimensionality, 3D crystal structure or sys- 
tem size. Our QMC results thus show that this asymp- 
totic \asge-U scenario is indeed realized in the finite-?/ 
parameter range relevant to the cuprates. We have also 
shown that on-site overscreening implies singularities in 
the imaginary-frequency dependence of the irreducible 
polarization insertion and of its underlying 3-point ver- 
tex function. 

When analyzed as an effective pairing potential i.e. 
formally, as a 1st order approximation to the irreducible 
particle-particle vertex, the screened potential Vs gives 
rise to a 1st neighbor d^^-y^ pairing attraction in the 
doping range between x > 5% and a; < 20 — 25% and to 
an on-site s-wave pairing attraction for x > 15%. The 
doping dependence of the d,j.2_y2 Eliashberg A-parameter 
closely tracks that of the extended (1st neighbor) attrac- 
tion of Vs, with a maximum of pairing strength near 
a; = 10— 15%, reminiscent of the observed doping depen- 
dence of the superconducting Tc in the cuprates. The 
doping dependence of the on-site s-wave A-parameter 
closely tracks that of the on-site Vs, increasing mono- 
tonically with x and becoming strongly attractive for 
X ~ 15%. This suggests the possibility of a doping in- 
duced transition or cross-over from d- to s-wave pairing. 
We have argued that retardation plays an essential role in 
the on-site overscreening of the interaction potential and 
that this should be refiected in the temporal / frequency 
structure of the corresponding s-wave superconducting 
order parameter. 
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